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ABSTRACT 

The renormahzation group flow in a general renormalizable gauge theory with a simple 
gauge group in 3+1 dimensions is analyzed. The flow of the ratios of the Yukawa couplings 
and the gauge coupling is described in terms of a bounded potential, which makes it 
possible to draw a number of non-trivial conclusions concerning the asymptotic structure 
of the theory. A classiflcation of possible flow patterns is given. 
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1 Beta functions of a general gauge theory 

When mass effects can be neglected, the behaviour of a quantum field theory under scale 
transformations can be described conveniently in terms of "running couplings" defined 
via beta functions. Whereas it is quite straightforward to calculate the beta functions, it 
is usually not so simple to deduce from them the structure of the renormalization group 
flow, especially when there is a large number of coupling constants. In the present work we 
will show that even for a very general theory it is possible to make a number of nontrivial 
statements concerning its asymptotic behaviour. 

We consider a gauge theory with a simple gauge group G and with Weyl fermions ipi 
and real scalars (pa transforming under some (generically reducible) representations Rp 
and Rb of G. The Lagrangian of this model is 



(1) 

where we have omitted gauge fixing, ghost, cjj^ and mass terms. Gauge invariance forces 
the Yukawa couplings Y^ and the quartic scalar couplings V"'^^^ to fulfill the relations 

Y,%Tf)1, + Y,tiTp)l^ + YSiTBf^^ = (2) 

and 

respectively. Tp and Tb are the generators of Rp and Rb- These equations could be 
solved by decomposing the group representations according to ipi — ^ ip-. , 0a — ^ 0^, where 
J, A run through the sets of irreducible representations spanned by z = z(J) and a = a{A). 
Then the couplings can be written as 

Y^^YAi]= T.iZ^'^)fji^'% (4) 

k 

where the A'^'^^ are tensors satisfying an analogue of Eq. (Q) whereas the {Z^''^)fj are the 
usual unrestricted coupling constants. In the same way also the \/°^'>'^ can be decomposed. 
It will turn out, however, that for a general analysis it is more convenient not to make 
use of this decomposition. 

The one- and two-loop contributions to the beta functions of this model have been 
calculated in the R^ gauge in dimensional regularization with the MS scheme |l[]. The 
gauge beta function is given by 

[QQf + QQb + CgilQSp + Sb- Mcg)] + 0{g'). (5) 



3(47r) 



C = T°'T°' is the quadratic Casimir Operator, dg the dimension oi G, S = j-Ti C the 
Dynkin index, and Q = j-TiC^. The counting of orders in g is such that Y = 0{g), 
V = 0{g^). The other beta functions are 

+r^Tr(r"^F^ + r'^tya) _ 6^2^^°^'^ + C^F")] + 0{g^ 
(using matrix notation for the fermionic indices) and 



Xs 



C4^\2^/'ya/37<5\ _ yapXeyySXe _|_ yaiXeyPSXe _|_ yaSXey/S-y 

+ [- Tt(Y"^Y^ + r^V") - Sg^C^^jV^/^^^ 
+ [- Tr(F^V^ + r^V^) - 3^2^^^]v-"A75 

+ [- Tr(F^^r^ + r^V^) - 3g^Cl^]V''f^^^ 

+ [- Tr(r^V^ + y^V^) - 3/C^^]l^"^^^ 

+3/({T^, T|r^{T^, T|r^ + {T^, TlJ^^T^, ^D^^^ + {T^, T|r^{T^, T|}^^) 
—2 Tr\(Y°''^Y^ + y/^^y")(y"^ty5 _|_ y<5ty7\ 

^(y«ty5 _^ y5ty.a^)^y/3ty7 + y7ty/3)] _^ c'(/). (7) 

In order to solve the system of ordinary differential equations (|^) - (|^) it is convenient 
to change variables to t = ln(/////o), G = [g/in)'^, y = Y/g and v = V/g'^. Then the 
evolution of G is governed by 



dG 2G^ I 3 

-^ = -AG2 - -^ 34c2 - lOcgSp - CgSB - QQf '^Qb + j 



2 _ ^G_ / 2 _ in^ c„ _ ^ c„ _ fin„ _fin„ -L A Tri'r'^.,/3^,,/3n _l di'/^^n 



where 

A = (22c, - 4^^ - Sb)/3. (9) 

To lowest order the solution is given by 

G-\t) = GQ^ + Xt + 0{G), (10) 

where Gq = G{t = 0) (similar notation will be used again), y and v evolve according to 

^ = |(4yW+yVV^+2/VV"+/Tr(y"t/+/t^.)_^«^_^^.)_^Q(^2)^ (11) 
dt 2 

where 

A, = 6(C^),, - A(5,,/2, (12) 



and 

(it 
+ [1 Tr(y^t^A ^ ^At^7) _ 3C^A ^ A^7A]^./3A. 

+ {y"^y' + /V")(/V^ + l/'^V)] + 0(G), (13) 

respectively. The forms of these equations suggest to define a new evolution parameter r 
by 

G 

dr = — dt and tq = 0. (14) 

r = -^ln(G'/Go) + 0(GlnG'), (15) 

2A 

r=^ + OiG'). (16) 



If A 7^ 0, then 



whereas for A = 



2 Potentials 

In the center of our further discussion of the flow of y there will be a "potential" U deflned 
by 

u = \t. \s^,>^l\" + E^"")' + E m,\' (17) 

ijkl 0/3 ij 

with 

s^Jkl = y%yti + yf.i/j'/ + ytiy%, (is) 

/"'3 = ImTri/°V, (19) 

and M = y''^y''-D. (20) 

It can easily be checked that, at lowest order in G, 

2 dyf, _ dyf^ _ dU 



G dt dr dy: 



a t ' 



(21) 



implying 



dU dU dyf^ , dU dyf^ ^,dU ^ 



^ = 21- — M>0 (22) 



dr (9|/"-^ dr dyf. dr dyf^ 



Thus f/ increases or decreases with increasing or decreasing r, along paths of steepest 
ascent or descent, respectively. U plays a role similar to that of c in Zamolodchikov's 
c-theorem p|. In fact, potentials for the renormalization group flow have been considered 
already 20 years ago 0. In general one has to define a metric with respect to which 
a gradient flow is defined. In our particular case it turns out to be just the Euclidean 
metric on y-space. The sign of U is chosen in such a way that /i, t, r, U are monotonically 
increasing functions of one another. 

From the explicit form of Eq. (|T7p it is obvious that U is non-negative, and it is 
not hard to show that the entries of M become large when one of the y^ becomes large. 
Therefore U has a global minimum, i.e. an infrared fixed point, at some finite value of y. If 
Rp or Rb contains isomorphic irreducible representations more than once, Eqs. (p^lf ) and 
(|13D are invariant under some global symmetries relating these irreducible representations 
and the set of fixed points generically will not be a unique point, but rather the orbit of 
such a symmetry. In fact it is also possible that the space of minima of U has a degeneracy 
that is not related to a global symmetry. An example is given by the model of Ref. [Q, 
where a continuous set of solutions of [/ = was found. Since the eigenvalues of Ti y°'^y^ 
would be invariant under global symmetries, solutions with different eigenvalues that were 
found cannot be related by symmetries. It would be interesting to find out whether local 
minima at different values of U are possible and whether the space of global minima is 
connected. 

It will also be important how fast the minimum is approached in the limit r ^ oo. 
If the Hessian if of f/ at ypp has maximal rank (i.e. if ypp is a non-degenerate critical 
point), all components of Ay = y — ypp will tend to zero exponentially in r because of 

^ = HAy + 0{{Ayn (23) 

As we have seen, the minimum often has flat directions, implying that the Hessian cannot 
have maximal rank. If the minimum locus is locally a submanifold of y-space (this is not 
fulfilled, for example, for the origin for U = ylyl)-, and if the rank of the Hessian is the 
maximal rank N (the dimension of y-space) minus the number N^a.t of fiat directions, 
one can locally define coordinates by the following procedure: Define A^fiat coordinates 
parametrizing the minimum locus and label each point in some suitably chosen open 
set by the A^'flat coordinates of the fixed point to which it flows and by A^ — Afiat linear 
combinations of the original coordinates on which the Hessian is positive definite. Under 
the flow, the first Afjat coordinates do not change, whereas the other A^ — Aflat coordinates 
approach the fixed point exponentially in r. Thus ||Ay|| is again exponentially bounded. 
If Aflat +rank(if) < A^? we cannot apply these arguments. For example, if dAy/dr ~ 
const. X (Ay)^, then Ay ~ const, x [const. — r)~^/^. 

We have written U in terms of three highly symmetric structures which were discovered 



in Ref. (see also 0) in the context of a search for finite theories. A case of special 
interest occurs when each of these structures vanishes separately. In particular this is the 
case for one loop finite supersymmetric theories ||^, 0]. For any supersymmetric theory 
Sijki will vanish: The Yukawa couplings come in two types, namely interactions involving 
gauginos, which are proportional to generators of the gauge group, and interactions within 
the chiral sector, which are determined by totally symmetric constants dijj^. Within one 
type of couplings contributions corresponding to the real and imaginary parts of the scalars 
cancel, whereas the mixed terms vanish because of the invariance condition on dijk (their 
contribution to Sijki is just the analogue of the l.h.s. of Eq. (^). The conditions / = 
and M = are not fulfilled automatically and lead to the well-known one-loop finiteness 
conditions for supersymmetric theories |P . An example for a non-supersymmetric theory 
with simultaneous vanishing of S^jki, 1°'^ and Mij is given in Ref. Q. There are many 
models, however, where it is not possible to solve f/ = 0. 

It is easy to see that U has at most one local maximum: Consider the one parameter 
family of couplings y"j{z) = z{yo)"j with a fixed set of values {yo)ij. Then U{y'^j{z)) is of 
the form a|z|^ + 6|2;p + c with a > 0, which has at most one local maximum, namely at 
z = 0. Therefore any maximum of U can only be located at yfj = for all a, i,j. Clearly 
y = is always a critical point. The Hessian there is diagonal (in the ds x (i|n-dimensional 
y-space) with diagonal elements corresponding to ?/j" given by —{Da + Djj). This point 
really represents a maximum if and only if the sub-matrix of D corresponding to Yukawa- 
interacting fermions is positive definite. If this is the case, y = is approached according 
to 

y'^j ~ const. X e-(^"+^^^)^ (24) 

(no summation over repeated indices). 

In addition to this maximum and the set of minima, U might (and usually will) also 
have saddle points. Particular examples are points where the Yukawa couplings for some of 
the particles vanish while the Yukawa couplings for the other particles correspond to some 
minimum of the restricted U . Although saddle points are unstable fixed points of the flow 
both in the infrared and in the ultraviolet limit, getting close to such a point (and thereby 
getting small beta functions) can considerably decelerate the flow, such that a realistic 
theory will often be close to a saddle point along large portions of the flow. A theory 
flowing exactly into a saddle point would require exactly fine-tuned initial conditions. 
In the absence of reasons for such a fine-tuning, this will occur with zero probability. 
A possible reason might be a global symmetry which could impose constraints on the 
couplings. In such a case one should consider the potential U only over the subspace of y- 
space defined by these constraints. Then the couplings run to extrema of the constrained 
potential. A typical example is supersymmetry which prevents y from running to in the 
ultraviolet limit. 



Let us now discuss the behavior of v: We consider Eq. ([T3|) for fixed y = ypp, neglecting 
higher orders in G. It is possible to integrate to get a potential again: 

U^ = v'"^^^v^^^'v^''''^ + 2[- Tr(y°V^ + y^V°) - 3C^^ + ^S''^]v'^'^^^v^P^^ 



+ (y"t/ + /V")(/V^ + y^V'^))]^;"^^' (25) 

In contrast to the potential for the Yukawa couphngs, U^ is unbounded, because it is 
cubic. We can show, however, that for any given set of y's it has at most a single local 
minimum: Assume there are two minima fi,W2 and consider the line v{s) = f i + sAf 
with Af = f2 — fi- Then s = and s = 1 are minima of the function f/„(f (s)) which 
is (at most) cubic in s. Therefore U^ must be constant along v{s). Now consider v{s) 
near s = 0. If every neighborhood of s = contains points that are not minima of (the 
full) f/„, then every neighborhood of Vi (in the full f-space) will contain points v' with 
U{v') < U{vi), in contradiction to the assumption that vi is a local minimum of Uy. 
If, on the other hand, there is a neighborhood of s = containing only minima of Uy, 
then d{vi + sAv)/dT must vanish identically. The coefficient of s^ in d{vi + sAv)/dT is 
proportional to At;"^-^^At;^^^^ + At;"^^'^At;^'^^^ + At;"^^^At;^'^^^ Summation over a = 7 
and (3 = 5 gives 

2 Y. {^V'^^^f + Y.CE At;""^-^)' = 0, (26) 

0/375 7(5 a 

implying At> = 0, i.e. vi = V2. The same arguments can be used to show that there is at 
most a single local maximum. 

We are now able to start a detailed case by case discussion of the asymptotic behavior 
(in the regime of validity of perturbation theory) of a theory described by the Lagrangian 
(|l]) both in the infrared and in the ultraviolet limit. 



3 Classification of fiow patterns 

A) A>0 

When r — >■ —00, G becomes large and perturbation theory is no longer valid. So, if there 
is any fixed point within the range of validity of perturbation theory, it must occur for 
r -^ +00, where G ~ e"^^"^. If the y's escape attraction by the ultraviolet fixed point or if 
there is no ultraviolet fixed point, then the y's go to infinity. One might wonder whether 
Y = gy could still remain finite. That this is not so follows from the fact that the y's 
reach infinity within finite r: 

— ~ const. X y^ implies — y^^ ~ const, x (r — const.) (27) 

Thus perturbation theory breaks down. 

If the starting point is in the region of attraction of the local maximum of U, all yg 
will go to zero. Near the fixed point, 

y^^ ~ const. X e-^^^^'^^"'^ ~ const, x g^d,'^^,,)/{2X) ^^8) 



Then, at lowest order in G, the evolution of v will be determined by Eq. (|1^) at y = 0. A 
simple example for a non-trivial fixed point of v is the case of just one real scalar which 
is a singlet under the gauge group: At y = 

G~^^ = 3v^ + Xv. (29) 

at 

There is an ultraviolet fixed point at fpp = —A/3 = — (22cc, — 4S'f)/9, which is approached 
according to 

V — vpp ~ const. X e^^^'^ ~ const, x G. (30) 

In general v — fpp will behave like some positive power of G. Describing the evolution of 
couplings in terms of some other coupling is just the idea of the "reduction of coupling 
constants" program |T^. The exponents of G encountered in y and v — fpp, which are 
not necessarily integral, correspond precisely (modulo the fact that different models were 



considered) to the non-integral exponents found in |]T0|. 



On the other hand there are many theories without a fixed point for v, due to the 



following argument (adapted from Ref. [TT|): A little calculation shows that at y = 






(31) 
which is positive definite for many theories. In such a case, or whenever we start outside 
the domain of attraction of an ultraviolet fixed point, for large r 

dv 

-— ~ const. X f^, v~^ ~ const, x (const. — r), (32) 

dr 

i.e. we leave the region of validity of perturbation theory within finite r. 

B) A<0 

Here G becomes large for r — *• +cxo, so that a perturbative fixed point is possible only in 
the infrared limit r — *■ — cxd. In this case y will always approach some infrared fixed point. 
According to the discussion after Eq. (^Bf), if A^flat+rank(if) is equal to the dimension of 
y-space, y — y-pp behaves like some exponential of r and therefore like a positive power of 
G. If this is not the case, it is possible that 

1/ - Z/FP ~ const. X [const. - t^^^^ ~ const, x (ln(G'"-^))"-^/^, (33) 

which is certainly a very slow approach to the fixed point. If a fixed point for v (at 
y = 2/Fp) exists and the starting values are in its domain of attraction, v will also go to 
the fixed point. 

As an example consider the case of A^ chiral fermions in some complex representation R 
of the gauge group and as many fermions in the conjugate representation R, together with 
a single scalar singlet. By a biunitary transformation one can diagonalize the couplings. 



so that i/ij = Zi(^i)5fj (cf. Eq. (|)), where I runs from 1 to A^. Then U turns out to be of 
the form 

U = aiY: \Ziff + hY.\Zi\'- cY: \Zi\^ + d (34) 

/ / / 

with positive constants a, 6, c, d. U is easily minimized with the resuh that at the fixed 
point all |2'/p are equal to some constant depending on A^ and the dimension and Casimir 
eigenvalue of R. At y = ypp, dv/dr is positive for v —>■ ±oo and negative for v = 0, i.e. it 
must have two zeroes. The larger of the two values of v for which dv/dr is zero corresponds 
to a minimum of U^. 

C) A = 

The behaviour of G is dictated by 



g^2 ^^ 2c^ + c.^F + Qf + Qb-^ TTiCpy^'y^) + 0(G), (35) 

where we have used A = 0, i.e. Sb = 22cg — 4S'i?, in order to eliminate Sb- y evolves 
according to Eq. (0) with D = GCp- 

Let us first consider the ultraviolet limit r ^ +oo: If no Yukawa interacting fermionic 
singlets are present, ?/ = is a local maximum of U. If there are fermionic singlets, or if 
the starting configuration is outside the domain of attraction, y will go to infinity within 
finite r. Let us now assume that perturbation theory remains valid for r —>■ +oo. Then 



liuiT-^+oo y = and stability of G, determined by Eq. (35) at y = 0, implies Qb < 2c^ 
and Sp < 2cg. Reinserting the latter inequality into A = 0, we get Sb > 14cg. Putting 
this into Eq. (|3l|) , we see that ti""TT will go to infinity within finite r. We conclude that 
for A = there is no perturbative ultraviolet fixed point. Of course, all this is again only 
true if there is no exact fine tuning of the initial values which would allow the couplings 
to stay in the minimum or in a saddle point of the potential. Such a case is considered in 
Ref. 0. 



For r -^ — oo, y will go to some infrared fixed point. The behavior of G is determined 
by Eq. ( pSf ) with Tr(Ci?|/^ y^) evaluated at the one-loop fixed point for y. If U has flat 
directions, ypp depends continually on the initial conditions. If all fiat directions corre- 
spond to symmetries of U, or if the fixed point is at U = 0, Tr^Gpy^ y^) is nevertheless 
stable under small changes in the initial conditions. If the r.h.s. of Eq. ( P5| ) is negative, 
perturbation theory will break down; if it is positive, G = will be an infrared fixed 
point, approached as G ~ {const. — t)^^/^. If U takes its minimum at vanishing Sijki, /"^ 
and Mij, then Tr(C^j/^V^) = eTrCj, = GdgQp and 

-2cl + CgSF-2QF + QB + i^.O. (36) 



8G^ dr 



The behavior of v depends on whether there is a fixed point for v at y = ypp and on the 
initial conditions. Of course the theory is stable in the infrared limit only if v approaches 
such a fixed point. 



The most interesting case occurs when the r.h.s. of Eq. (|35D, evaluated at y = ypp, 
vanishes. Then 

^,§ = MCiy^'y^ - .p^pVf^p)) + 0(G) 

= Tr(C^(A/^A/ + y^p^A/ + A/V^p)) + 0(G), (37) 

where Ay = y — ypp, so it is crucial for the further discussion how the fixed point is 
approached. If there are components of y behaving like (const. — r)^^/^, G will run to 
or oo. If, however, A^flat+rank(i7) is equal to the dimension of y-space, all components 
of Ay will tend to zero exponentially in r. Again the same type of behaviour has been 
encountered in the context of "reduction of coupling constants" [|13]. The r.h.s. of Eq. 
( |37D receives an exponential inhomogeneity; for suitable (but not exactly fine-tuned) initial 
values G remains finite for r — > — oo at the present order in perturbation theory. In a 
next step one has to consider higher order corrections to the evolution of y. Using the 
vanishing of the first two orders oi (3{g), one gets 

with /3(2)(r) = gG'^f{y, v). The ansatz y = ypp + Gy'^^^ + ... yields 

G^ = G{Hy^^^ + 2/(2/FP, t^Fp)) + h.O. (39) 

Thus y^^^ goes to —2H^^f{ypp, vpp) exponentially with exactly the same rate as Ay goes 
to zero. In a similar way fpp gets shifted by a term of the order of G. The value of y^^^ 
is relevant for the G^-term (which is then dominant) in f3{G). The same analysis can be 
repeated until eventually G starts to run like some power of r. One would expect that 
this should always happen at some order in perturbation theory. This is however not the 
case for supersymmetric theories, due to the theorem of Ref. |T^ that in such theories 
the vanishing of all beta functions at A^ loop level implies vanishing of the gauge beta 
function at A^ + 1 loop level. 

Summing up our results, we can distinguish the following four cases: 

• Perturbation theory remains valid neither in the infrared nor in the ultraviolet limit. 

• The theory is unstable in the infrared limit but shows asymptotic freedom: For 
t -^ oo, g behaves according to 5^ ~ t~^^^. The Yukawa couplings go to zero 
as powers of g with exponents greater than 1; the behavior of the quartic scalar 
couplings is determined by a nontrivial fixed point of V/g"^. 

• The theory is unstable in the ultraviolet limit. In the infrared limit g goes to zero, 
Y/g and V/g'^ approach nontrivial fixed points. 

• There is no ultraviolet stable fixed point; in the infrared limit all couplings approach 
a fixed point. 



There is no theory with stable fixed points both in the infrared and in the ultraviolet 
limit. 
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